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EIGENMODE ANALYSIS OF BOUNDARY CONDITIONS 
FOR THE ONE-DIMENSIONAL PRECONDITIONED EULER EQUATIONS 

DAVID L. DARMOFAL* 


Abstract. An analysis of the effect of local preconditioning on boundary conditions for the subsonic, 
one- dimensional Euler equations is presented. Decay rates for the cigenmodcs of the initial boundary value 
problem are determined for different boundary conditions. Riemann invariant boundary conditions based 
on the unprcconditioncd Euler equations are shown to be reflective with preconditioning, and, at low Mach 
numbers, disturbances do not decay. Other boundary conditions arc investigated which are non- reflective 
with preconditioning and numerical results are presented confirming the analysis. 

Subject classification. Applied & Numerical Mathematics 
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1. Introduction. Local preconditioning has been successfully utilized to accelerate the convergence to 
a steady-state for Euler and Navicr-Stokes simulations [15, 18, 2, 19, 17, 7, 11, 4, 3]. Local preconditioning 
is introduced into a time-dependent problem as, 

u t + P(u)r(u) = 0, 


where u is the state vector of length m, r is the residual vector of length m, and P is the m x m pre- 
conditioning matrix. Since preconditioning effectively alters the time-dependent properties of the governing 
partial differential equation, modifications of the numerical discretization can be required. For example, 
upwind methods for inviscid problems must be based on the characteristics of the preconditioned equations 
instead of the unpreconditioned equations[18]. Similarly, the behavior of boundary conditions in conjunction 
with preconditioning will also be altered. While the effect of preconditioning on boundary conditions is 
known[9, 6, 16], to date, no quantitative analysis has been performed. 

The purpose of this paper is to analyze the effect of preconditioning on several different boundary condi- 
tions commonly used in numerical simulations. Specifically, we consider the one-dimensional, preconditioned 
Euler equations linearized about a steady, uniform, subsonic mean state. The work is an extension of the 
analysis of Giles [5] for the one-dimensional, unpreconditioned Euler equations. As discussed by Giles, the 
exact eigenmodes and cigcnfrequencies for this initial boundary value problem can be analytically deter- 
mined. From these, we find the exponential decay rates for perturbations under different sets of boundary 
conditions. Finally, we demonstrate the validity of the analysis through numerical results. 

2. Analysis. We first present the general analysis of the initial boundary value problem following the 
work of Giles[5]. The linear, preconditioned Euler equations are given by, 


( 2 . 1 ) 
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where p , q , p are the perturbation density, velocity, and pressure, and p, q , c arc the mean density, velocity, 
and speed of sound. The speed of sound is related to the pre ssure and density through c 2 = y Pi P* The 
subsonic inflow is located at X — 0 and the outflow is at X — h. 

Next, we define the following non-dimensionalizations to simplify the analysis, 


( 2 . 2 ) 



The non-dimensional version of Equation (2.1) is, 

(2-3) Ut + PAuj = 0, 


where, 
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and the mean flow Mach number is M = q/c. The boundary conditions for subsonic flow require two inflow 
quantities and one outflow quantity to be specified. The inflow boundary conditions can be expressed as, 


(2.4) Cj n u(0, t) = 0, 

where C m is a 2 x 3 matrix dependent on the specific choice of inflow conditions. Similarly, the single outflow 
boundary condition can be expressed as, 


(2-5) C owt u(M) = 0, 

where C ou t is a 1 x 3 matrix dependent on the specific choice of outflow conditions. 

Equations (2.3), (2.4), and (2.5) represent the initial boundary value problem which we wish to study. 
An eigenmode of the initial boundary value problem is given by 

3 

(2.6) u = e ~ iu;t <*iTie lU)x/Xi . 

i= 1 

r i and A* are the right eigenvectors and eigenvalues, respectively, of the matrix PA, i.e., 

(PA - AJ) r* = 0. 

In the following developments, we assume that the eigenvalues have been ordered such that the two forward- 
moving characteristics arc i = 1,2 (i.e. A 12 > 0) and the backward- moving characteristic is i = 3 (i.e. 
A 3 < 0). The constants a t are the strengths of each eigenmode. Given an eigenmode of the initial boundary 
value problem as in Equation (2.6), the strength can be determined by multiplying u by the left eigenvector, 
1, , where left and right eigenvectors are related by, 

(2-7) ( li b 1 3 ) = ( ri r 2 r 3 j 

Specifically, since Equation (2.7) implies that \frj = Sij, we may find, 

lfu = 1 Te~ lu;t J2 a j r j eiuJx/Xj = o 
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The eigenfrequency u) and eigenmodc strengths a* arc determined by the boundary conditions. For the 
inflow boundary, substitution of Equation (2.6) into Equation (2.4) leads to, 


( 2 . 8 ) 


where 


(2.9) 
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As described by Giles[5], a necessary condition for the well-posedness of the initial boundary value problem is 
that the incoming characteristics, ot\ and can be determined as functions of the outgoing characteristic, 
a 3 . This requires that the 2x2 matrix, 



is non-singular. Also, the boundary condition at the inflow will be non-reflecting if the outgoing characteristic 
does not cause a perturbation in the incoming characteristics. In other words, a necessary condition for non- 
reflecting inflow conditions is that &i 3 = 623 = 0. 

For the outflow boundary, substitution of Equation (2.6) into Equation (2.5) leads to, 


( 2 . 10 ) 

where 
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In this case, well-posedness of the initial boundary value problem requires that incoming characteristic, a 3 , 
can be determined as a function of the outgoing characteristics, oc\ and a^* Thus, 6 33 must be non-zero. 
Also, the boundary condition at the outflow will be non-reflecting if 6 3 \ = 6 3 2 = 0. 

The inflow and outflow boundary conditions in Equations (2.8) and (2.10) can be combined as, 


( 2 . 12 ) 


B(w 



0. 


In order for a non-trivial solution of the initial boundary value problem to exist, a non-zero vector, (aj , <22, &3) T , 
must exist which satisfies Equation (2.12). This is possible for values of lo for which, 


det B(cj) = 0. 


Separating the eigenfrequency into its real and imaginary parts, uj = cj r + to;*, the amplitude of the eigenmodes 
grows as exp(cjit). Thus, in order for the eigenmodes to decay, < 0 for all possible values of lj. 

Finally, we note that the steady-state problem is well-posed if and only if detB(0) is non-zero. [5]. 
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3. Examples. In this analysis, wc will consider the one-c imcnsional version of the van Leer-Lee-Roc 
preconditioner [ 1 8 ] which is discussed by Lee[ 8 ]. For this preconditioncr, the resultant PA matrix is, 
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The specific form of P is described in the Appendix. The eigenvalues of PA are Ai , 2i3 = M, M, -M, and 
the eigenvectors are, 


( 3 . 1 ) 
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Some physical interpretation can be given to the preconditioned eigenmodes by considering the strengths 
of each mode in dimensional form, i.e., 


( 3 . 2 ) 
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( 3 . 3 ) 
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The strength of the first eigenmode is proportional to the linearized change in entropy and, under this 
preconditioning, is unchanged from the Euler equations without preconditioning. The second eigenmode 
is proportional to the linearized streamwise momentum, p + paq. Finally, the third eigenmode is directly 
proportional to the pressure. Thus, with preconditioning, the upstream- running wave is exactly a pressure 
wave. As we will show in Section 3 . 2 , this will have a significant impact on the behavior of boundary 
conditions for which a downstream pressure is specified. 

3.1. Euler Riemann boundary conditions. We first consider the use of boundary conditions based 
on the Riemann invariants of the Euler equations. Without preconditioning, Giles[ 5 ] showed that these 
boundary conditions are non-reflecting with w* = - oo. In other words, all perturbations are eliminated 
in the finite time required for them to propagate out of the domain. Specifically, the Riemann invariant 
boundary conditions arc, 


A = 0, 


/ P'/P n 

l * + 
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where the primed quantities are the sum of the mean and pertu bation states, i.e. p' = p + p. Linearization 
and non-dimensionalization of these boundary conditions gives, 
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C out = ( 1 7-1 -7 j . 
/.From Ci n and C out , we find the matrix B, 
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The determinant of B is, 


det B = —(7 - l) 2 


(M - + {M+ l)e“ iw/M 


and the eigenfrequencies which result in a zero determinant arc, 


uj r = nMn , for integer n, 
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In contrast to the unpreconditioned Euler equations, the Riemann invariant boundary conditions are reflective 
for the preconditioned Euler equations since the value of Ui is not —00. Within a computational simulation, 
perhaps the best measure of the decay rate is actually the decay per time step. Assuming the time step is 
given by a CFL condition of the form A t = ^Ax/A max where v is a constant dependent on the temporal 
integration, Ax is the cell size, and A max is the maximum eigenvalue, then u)At oc ^/A max . Since A max = M 
for this preconditioned system, 


1 x _ It 1 + M 

Amax — 2 1 — M 

In particular, we note that as M — > 0, u>i/ A max — > 0. Thus, at low Mach numbers, disturbances will not 
decay indicating that the use of Riemann invariant boundary conditions based on the Euler equations is 
likely to impede convergence to a steady state. 


3.2. Entropy, stagnation enthalpy at inflow; pressure at outflow. Another common set of 
boundary conditions for subsonic, internal flows is the specification of entropy and stagnation enthalpy at 
the inflow and pressure at the outflow. Specifically, 


A = 0, 
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For these boundary conditions, C i n and C ou t are, 

' 0 1 
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^From C tn and C 01i £, we find the matrix B, 


-10 0 \ 

-1 ( 7-1 )M 0 

0 0 Me~ iu/M ) 


We note that 613 = 623 = ^31 = 632 — 0 for these boundary conditions. Thus, outgoing waves do not 
generate incoming waves at either the inflow or outflow boundaries. At the inflow, specification of the 
entropy guarantees that 612 = 613 = 0 since the change in entropy is proportional to the strength of the 
first eigenmode as shown in Equation (3.2). Specification of the stagnation enthalpy is also a non-reflecting 
inflow condition since the change in stagnation enthalpy, H , can be written as a linear combination of the 
first and second cigcnmodes, H oc — li T u + (7 — 1)M1 2 T u. Finally, as observed from Equation (3.4), the 
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pressure change is proportional to the change in the third (upsti earn- running) cigenmode, thus, specification 
of the pressure at the outflow is non- reflective. 

Since the boundary conditions at inflow and outflow arc individually non- reflective, the entire system 
will have infinite, perturbation decay rates. Specifically, the de erminant of B is, 

det B = -(7- 1 

and the eigenfrequcncics which result in a zero determinant are, 


<jj r = 27 rMn, for integer n, 

<jJ x — —00. 


Since uj x — —00, disturbances are eliminated via propagation In finite time. An interesting aspect of this 
result is that these boundary conditions are actually reflective for the unpreconditioned Euler equations (sec 
Giles[ 5 ]). 

3*3. Velocity, temperature at inflow; pressure at outflow. The final set of boundary conditions 
we consider is setting velocity and temperature at the inflow arid pressure at the outflow. These conditions 
are fairly common in low speed viscous flow applications. Specifically, 


P f /P f = p/p 


q = q 


X = 0 , 

For these boundary conditions, C m and C ou t arc, 

C. n = ^ 

Ccmi = ^ 

^From C m and C outl we find the matrix B, 
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Unlike the boundary conditions in Section 3 . 2 , the specification of velocity and temperature at the inflow is 
a reflective condition and create reflections of the outgoing chara :teristic wave. This is evident from the non- 
zero values of 613 and 623. Specifically, at the inflow, oc\ = (7 — l)Mcx^ and 0,2 = 0:3. However, since 631 = 
^32 = 0 . the outflow boundary condition does not create any perturbations in the incoming characteristics. 
Thus, the reflected waves from the inflow conditions would be emitted at the outflow boundary without 
further reflection and all disturbances would be eliminated in finite time. Specifically, the determinant of B 
is, 



det B = 

and the eigenfrequencies which result in a zero determinant are, 

bj T — 2 k M n, for integer n, 

U)i = -00. 
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While perturbations again have infinite decay rates for these boundary conditions, in practice, the reflec- 
tive inflow condition may slow convergence somewhat compared to the non-reflecting boundary conditions 
discussed in Section 3.2. This convergence slowdown is observed in the numerical results of the following 
section. 

4. Numerical Results. To illustrate the effect of different boundary conditions on numerical conver- 
gence as well as check the accuracy of the analysis, we simulate the two-dimensional flow in a duct with a 
straight upper wall and a bump on the lower wall between 0 < x < 1 described by y = .042 sin 2 (7r:c). The 
domain is 5 unit lengths long and 2 lengths high. The grid is structured with clustering toward the wall 
boundary. 

Specifically, we use the algorithm described by Darmofal and Siu[3] which employs the semi-coarsening 
technique of Mulder [12, 13] in conjunction with a multi-stage, block Jacobi relaxation[10, 1]. The dis- 
cretization is a 2nd order upwind scheme with a Roc approximate Riemann solver [14]. The calculations are 
performed on a grid of 32 x 16 cells. A three level, V-cyclc is utilized with 2 pre and post smooths. All 
calculations arc initialized to uniform flow. 

A form of Turkel’s preconditioner is employed which is smoothly turned off for Mach numbers above 0.5. 
While this preconditioner is different than the one- dimensional van Leer-Lee-Roe preconditioncr for which 
the analysis was performed, we expect the low Mach number behavior of these preconditioners to be similar. 
The major differences between the analysis and the numerical results will occur for higher subsonic Mach 
numbers where the preconditioning is turned off in the numerical simulations. 

We have implemented the boundary conditions described above by constructing a boundary face state 
vector and calculating the boundary flux directly from this state vector. For example, for the enthalpy, 
entropy, pressure boundary conditions at an inflow, entropy, enthalpy, and the tangential velocity are pre- 
scribed from the exterior and the pressure is extrapolated from the interior. At an outflow, we reverse the 
procedure and specify pressure from the exterior and extrapolate entropy, enthalpy, and tangential velocity 
from the interior. Note, regardless of the specific boundary conditions, we always use the tangential velocity 
as the additional variable for the two-dimensional boundary implementation. 

The number of cycles required to converge the solution six orders of magnitude from the initial residual 
are given in Table 4.1. As can be clearly seen, the results are in good agreement with the analysis at 
low Mach numbers. In particular, the Riemann boundary conditions are unstable at low Mach numbers. 
The entropy, enthalpy, pressure (HSP) boundary conditions perform best while the velocity, temperature, 
pressure boundary (QTP) conditions are about 75% more expensive. This would indicate that the reflective 
nature of the inflow for the QTP boundary conditions does indeed slow the convergence. 

At the higher Mach numbers (M > 0.2), the Riemann boundary condition cases begin to converge and 
the number of cycles decreases with increasing Mach. The QTP boundary conditions require an increasing 
number of cycles for increasing Mach number. Also, for the M 0 c =0.5 test case, the amount of work required 
to converge the HSP boundary conditions increased from 8 cycles to 10 cycles. These trends at higher Mach 
numbers arc expected because the preconditioning which was implemented numerically was automatically 
phased out at M = 0.5. Thus, the behavior of the different boundary conditions will be described by the 
unpreconditioned Euler equations. In this case, the Riemann boundary conditions are non- reflective while 
the other boundary conditions are reflective. 

5. Remarks. The present analysis of the preconditioned Euler equations shows the effect of precondi- 
tioning on boundary conditions. Boundary conditions based on the Riemann invariants of the Euler equations 
are found to be reflective in conjunction with preconditioning. The problem is most detrimental at low Mach 
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Mach 

Riemann 

HSP 

QTP 

0.001 

UNS 

8 

13 

0.01 

UNS 

8 

13 

0.1 

UNS 

8 

13 

0.2 

20 

8 

14 | 

0.3 

14 

8 

15 

0.4 

11 

8 

18 

0.5 

9 

10 20 


Table 4.1 

Number of cycles required to drop residual six orders of magnitude for different Mach numbers and boundary conditions. 
Riemann: Euler Riemann invariant boundary conditions from Section 3 1. HSP: enthalpy, entropy, pressure boundary con- 
ditions from Section 3.2. QTP: velocity, temperature, pressure boundary conditions from Section 3.3. UNS: algorithm was 
unstable and aborted with infinite residual. 


numbers where the decay rate of perturbations approaches zero. Boundary conditions which specify entropy 
and stagnation enthalpy at an inflow and pressure at an outflow are found to be non- reflective with precon- 
ditioning. Numerical results were presented which are in good agreement with the analytic predictions. 

Finally, an interesting possibility implied by this analysis would be to incorporate boundary condition 
considerations into the design of preconditioners. For example, [given a set of physical boundary conditions 
which must be applied for a specific problem, a preconditioner dould be designed such that these conditions 
are non- reflective. 


Appendix. The one- dimensional van Leer-Lee-Roe precon iitioner [8] is usually derived using the sym- 
metrizing variables which in dimensional form are (p/pc, q, p - c 2 p). Using Equation (2.2), the non- 
dimensional symmetrizing variables are, 


v = ( p, q, p- py , 

and arc related to the u = (p, q,p) T variables through the transformation, v = Su, where, 

/ 0 0 1 


s = 


0 1 0 
V -i o i 


The preconditioned Euler equations in terms of v arc, 

v« ~\~ P = 0, 

where, P v = SPS”\ and, 

( M 1 0 

1 M 0 
0 0 M 

Finally, the one- dimensional van Leer-Lee-Roe preconditioncr is given by, 
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and, for use with the (p,q,p) variables, 


where (3 2 = 1 — M 2 . 


P = 


f i M Af 2 , \ 

1 ~w w 1 \ 

0 i + £ -% 
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